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Abstract. We consider the effect on the propagation of light of inhomogeneities with 
sizes of order 10 Mpc or larger. The Universe is approximated through a variation 
of the Swiss-cheese model. The spherical inhomogeneities are void-like, with central 
underdensities surrounded by compensating overdense shells. We study the propagation 
of light in this background, assuming that the source and the observer occupy random 
positions, so that each beam travels through several inhomogeneities at random angles. 
The distribution of luminosity distances for sources with the same rcdshift is asymmetric, 
with a peak at a value larger than the average one. The width of the distribution 
and the location of the maximum increase with increasing redshift and length scale 
of the inhomogeneities. We compute the induced dispersion and bias on cosmological 
parameters derived from the supernova data. They are too small to explain the perceived 
acceleration without dark energy, even when the length scale of the inhomogeneities 
is comparable to the horizon distance. Moreover, the dispersion and bias induced by 
gravitational lensing at the scales of galaxies or clusters of galaxies are larger by at least 
an order of magnitude. 



Light Propagation and Large-Scale Inhomogeneities 



2 



1. Introduction 

The observed deviation from homogeneity of the structure of the Universe at small length 
scales poses the question of whether the use of the Friedmann-Robertson- Walker (FRW) 
metric is adequate for the discussion of the cosmological expansion. The argument for 
the applicability of a homogeneous solution is based on the observation that the matter 
distribution is homogeneous when averaged over length scales of 0(100) h~ l Mpc. On 
the other hand, the conclusion drawn from observations of distant supernova [1, 2], 
or the perturbations of the cosmic microwave background (CMB) [3], that the recent 
cosmological expansion is accelerating has placed the effect of structure formation on the 
overall expansion under scrutiny. 

We are interested in the influence of inhomogeneities with sub-horizon characteristic 
scales on the perceived expansion [4, 5, 6]. The observed inhomogeneities in the matter 
distribution with length scales of (9(10) h~ l Mpc or smaller are large, so that the Universe 
cannot be approximated as homogeneous at these scales. It is important to have a 
quantitative estimate of the influence of these inhomogeneities on the data used for the 
determination of the expansion rate. As all the observations involve the detection of light 
signals, it is crucial to understand the effect of the inhomogeneous background on light 
propagation. 

The transmission of a light beam in a general gravitational background can be studied 
through the Sachs optical equations [7]. These describe the expansion and shear of the 
beam along its null trajectory. Apart from the case of a FRW background, the optical 
equations have been derived by Kantowski for a Scharzschild background [8]. He used 
them for the study of light transmission within the Swiss-cheese model of the Universe [9]. 
In this model the light propagates essentially in empty space with a beam expansion larger 
than the average. In rare instances it passes near a very dense clump of matter, which 
produces significant shear and focusing of the beam [10]. Because of the randomness of 
such events, sources with the same intrinsic luminosity and redshift may have different 
luminosity distances. The distribution is peaked at a value of the luminosity distance 
larger than the one in a homogeneous background [11, 12], even though the mean value 
remains unaffected. If the sample of observed sources is small, it unlikely that a clump 
of matter will be encountered during the beam propagation. The beams essentially 
propagate in empty space, while the expansion is induced by the average energy density. 
A phenomenological equation has been proposed by Dyer and Roeder in order to describe 
the case in which a fraction of the energy density is homogeneously distributed while the 
remaining is in clumps that do not affect the light propagation [13]. 

The picture of the Universe we described above is applicable to length scales of 
(9(1) h~ x Mpc or smaller, for which the dominant structures are galaxies or clusters 
of galaxies. At larger distances the averaged matter distribution has a smaller density 
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contrast. The use of the Schwarzschild geometry for the description of the inhomogeneities 
is not appropriate. The Lemaitre-Tolman-Bondi (LTB) metric [14] has been employed 
often for the modelling of the Universe at scales of O(10) h~ x Mpc or larger [15]— [24] . Its 
use demonstrates that inhomogeneities can induce deviations of the luminosity distance 
from its value in a homogeneous background. For example, it has been observed that 
any form of the luminosity distance as a function of redshift can be reproduced with 
the LTB metric [15]. This means that the supernova data can be explained through an 
inhomogeneous matter distribution in the context of this metric. However, reproducing 
the data requires a variation of the density or the expansion rate over distances of 
(9(100) h~ l Mpc or even larger [20]-[24]. In order to avoid a conflict with the isotropy 
of the CMB the location of the observer must be within a distance of O(10) h~ l Mpc 
from the center of the spherical configuration described by the LTB metric. In this sense, 
the explanation of the perceived acceleration relies on the position of the observer. On 
the other hand, there are indications for the presence of a very large void in our vicinity 
[25, 26]. 

Our work is based on the fundamental assumption that we do not occupy a special 
position in the Universe. We are interested in determining the maximum effect that 
large-scale inhomogeneities can have on the luminosity distance. The density contrast 
of inhomogeneities with characteristics lengths of O(10) h~ l Mpc or larger is at most of 
(9(1). For this reason it seems unlikely that they can affect the propagation of light more 
strongly than the inhomogeneities with lengths of (9(1) h~ l Mpc or smaller, whose density 
contrast can be larger than 1 by several orders of magnitude. On the other hand, the 
validity of the FRW metric is questionable if inhomogeneities with lengths comparable 
to the horizon distance ~3x 10 3 /i _1 Mpc develop a density contrast of (9(1). We shall 
allow for such a possibility in order to obtain a quantitative estimate of the effect on the 
luminosity distance. 

In a previous publication [27] we derived the optical equations for a general LTB 
background [14]. We used them in order to study light propagation in a variant of 
the Swiss-cheese model. The inhomogeneities are modelled as spherical regions within 
which the geometry is described by the LTB metric. At the boundary of these regions 
the LTB metric is matched with the FRW metric that describes the evolution in the 
region between the inhomogeneities. The Universe consists of collapsing or expanding 
inhomogeneous regions, while a common scale factor exists that describes the expansion 
of the homogeneous intermediate regions. This model is similar to the standard Swiss- 
cheese model [9], with the replacement of the Schwarzschild metric with the LTB one. 
For this reason we refer to it as the LTB Swiss-cheese model. 

We focus on the effect of inhomogeneities on the luminosity distance if the source and 
the observer do not occupy special positions in the Universe. This is achieved by placing 
both the source and the observer within the homogeneous region of the LTB Swiss-cheese 
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model. During its path the light signal crosses several inhomogeneous regions before 
reaching the observer. As we have mentioned, similar studies [11, 12] have discussed 
the influence of structures such as galaxies or clusters of galaxies on light propagation. 
The beam shear plays an important role in this effect, characterized as gravitational 
lensing [28]. We focus on inhomogeneities of much larger length scale, of O(10) h" 1 
Mpc or larger. The matter distribution, even though inhomogeneous, is more evenly 
distributed than in the previous case. In particular, we assume that each spherical region 
has a central underdensity surrounded by an overdense shell. The densities in the two 
regions are comparable at early times and differ by a factor 0(1) during the later stages 
of the evolution. The LTB Swiss-cheese model we construct in this way describes a 
Universe dominated by spherical voids with the compensating matter concentrated in 
shells surrounding them. The beam shear is negligible in our calculations, and the main 
effect arises from the variations of the beam expansion because of the inhomogeneities. 

In ref. [27] we estimated the deviations of the luminosity distance from its value in a 
homogeneous Universe by considering the extreme case in which the light passes through 
the centers of all the inhomogeneities it encounters. In this work we perform a more 
detailed statistical analysis, by considering light beams with random impact parameters 
relative to the centers. We check that the resulting luminosity distance as a function of the 
impact parameter is consistent with flux conservation. For a given redshift we estimate the 
width of the distribution of luminosity distances. We also determine the deviation of the 
maximum of the distribution from the value in a homogeneous background. For a central 
underdensity this value is positive and sets the scale for the perceived increase in the 
expansion rate relative to a homogeneous Universe. An analytical study with similarities 
with our work is described in ref. [29]. It focuses, however, on the study of the redshift, 
which determines only partially the luminosity distance. No statistical analysis, which is 
the central point of our work, is performed either. 

In the following section we summarize the geodesic and optical equations in a LTB 
background. In section 3 we describe the cosmological evolution of the inhomogeneous 
regions. In section 4 we study, both numerically and analytically, the effect of the 
inhomogeneity on the properties (redshift, beam area) of a light beam that travels through 
it. In section 5 we study the luminosity distance as a function of the angle at which the 
beam crosses the inhomogeneity. We verify that the results are consistent with flux 
conservation. In section 6 we consider multiple crossings of inhomogeneities by the light 
beam and the effect on the luminosity distance. We calculate the distribution of the 
deviations of the luminosity distance from its value in a homogeneous background. In 
section 7 we estimate the effect on the determination of cosmological parameters from 
supernova data. 
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2. Luminosity distance in a Lemaitre-Tolman-Bondi (LTB) background 

Under the assumption of spherical symmetry, the most general metric for a pressureless, 
inhomogeneous fluid is the LTB metric [14]. It can be written in the form 

ds 2 = -dt 2 + b 2 (t, r)dr 2 + R 2 (t, r)dtt 2 , (2.1) 

where dfl 2 is the metric on a two-sphere. The function b(t, r) is given by 

*"M = g*$, (2,) 

where the prime denotes differentiation with respect to r, and f(r) is an arbitrary function. 
The bulk energy momentum tensor has the form 

Tl = diag (-/>(*, r), 0, 0, 0) . (2.3) 

The fluid consists of successive shells marked by r, whose local density p is time-dependent. 
The function R(t, r) describes the location of the shell marked by r at the time t. Through 
an appropriate rescaling it can be chosen to satisfy -R(0, r) = r. 
The Einstein equations reduce to 

M'(r) =47iR 2 pR', (2.5) 

where the dot denotes differentiation with respect to t, and G = (16ttM 2 ) . The 
generalized mass function A4(r) of the fluid can be chosen arbitrarily. Because of energy 
conservation Ai{r) is independent of t, while p and R depend on both t and r. 

Without loss of generality we consider geodesic null curves on the plane with 9 = 7r/2. 
The first geodesic equation is 

dk° R'R' / A 2 • / „x2 , s 

— + jik 1 ) + RR (k) =0, (2.6) 



dX 1 + f 

with k % = dx l /d\ and A an affine parameter along the null beam trajectory. The second 
geodesic equation can be replaced by the null condition 

-(A;°) 2 + I ^(fc 1 ) 2 + R 2 (k*) 2 = 0, (2.7) 

while the third one can be integrated to obtain 

k 3 = § ■ (2.8) 
The equation for the beam area can be written as [27] 

1 d 2 \J~A 1 /,n\ 2 9 / 

The shear a is generated by inhomogeneities, for which the local energy density is different 
from the average one [27]. It describes the deformations of the cross-section of the beam, 
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induced by the propagation within an inhomogeneous medium. The shear is important 
when the beam passes near regions in which the density exceeds the average one by several 
orders of magnitude. Within our modelling of large-scale structure, applicable for scales 
above O(10) h~ l Mpc, the average density contrast is not sufficiently large for the shear 
to become important. (We have verified this conclusion numerically.) For this reason we 
neglect it in our study. 

The Friedmann-Robertson- Walker (FRW) metric is a special case of the LTB metric 

with 

R(t,r) = a(t)r f(r) = cr 2 , c = 0, ±1 (2.10) 

P=^y M(r) = f V ». (2.11) 

The geodesic equation (2.6) has the solution 

k° = ^- y (2.12) 

The solution of eq. (2.9) for an outgoing beam is 

A(X) = r 2 (A) a 2 (t(X)) Q s . (2.13) 

If we normalize the scale factor so that a(t s ) = 1 at the time of the beam emission, 
we recover the standard expression A = r 2 Q s in flat space-time. The constant Q s can 
be identified with the solid angle spanned by a certain beam when the light is emitted 
by a point-like isotropic source. We are interested in light propagation in more general 
backgrounds. We assume that the light emission near the source is not affected by the 
large-scale geometry. By choosing an affine parameter that is locally A = t in the vicinity 
of the source, we can set 
aly/A 



This expression, along with 
C4 



n s . (2.14) 

A=0 



= 0, (2.15) 

A=0 ' V ' 

provide the initial conditions for the solution of eq. (2.9). 

In order to define the luminosity distance, we consider photons emitted within a 
solid angle fl s by an isotropic source with luminosity L. These photons are detected by 
an observer for whom the light beam has a cross-section A Q . The redshift factor is 

1+ * = S = f < 2 - 16 ' 

because the frequencies measured at the source and at the observation point are 
proportional to the values of k° at these points. The energy flux f a measured by the 
observer is 

l l n s 
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The above expression allows the determination of the luminosity distance Dl as a function 
of the redshift z. The beam area can be calculated by solving eq. (2.9), with initial 
conditions given by eqs. (2.14), (2.15), while the redshift is given by eq. (2.16). 

3. Modelling the inhomogeneities 

We study the effect of inhomogeneities on the luminosity distance without assuming a 
preferred location of the observer. We model the inhomogeneities as spherical regions 
within which the geometry is described by the LTB metric. At the boundary of these 
regions, the LTB metric is matched with the FRW metric that describes the expansion of 
the homogeneous intermediate regions. Our model is similar to the standard Swiss-cheese 
model [9], with the replacement of the Schwarzschild metric with the LTB one. For this 
reason we refer to it as the LTB Swiss-cheese model. 

The choice of the two arbitrary functions Ai(r) and f(r) in eq. (2.1) can lead to 
different physical situations. The mass function A4(r) is related to the initial matter 
distribution. The function fir) defines an effective curvature term in eq. (2.4). We 
work in a gauge in which i?(0, r) = r. We parametrize the initial energy density as 
pi{r) = (1 + e(r)) p 0>i , with Pi(r) = p(0,r) and \e(r)\ < 1. The initial energy density of 
the homogeneous background is po,i- If the size of the inhomogeneity is ro, a consistent 
solution requires 47r Jq° r 2 e{r)dr = 0, so that 



This is obvious if we apply eqs. (2.4), (2.5) to the homogeneous part at r > r . The FRW 
metric is a special case of the LTB metric, described by eqs. (2.10), (2.11). Consistency 
of the equations as the shell with r = r is approached from both sides imposes the 
condition (3.1). Moreover, the absence of singularities requires the continuity of f(r) and 
fir). As we assume that the homogeneous part is flat, this means that /(ro) = /'( r o) = 0. 
Alternatively, one may consider the matching of the two metrics at the surface with r = ro 
employing junction conditions [30]. If this surface does not contain a singular energy 
density, the above constraints must be imposed [29, 31]. 

We assume that at the initial time t{ = the expansion rate Hi = R/R = R'/R' is 
given for all r by the standard expression in homogeneous cosmology: Hf = p Qi /(6M 2 ). 
Then, eq. (2.4) with R(0,r) = r implies that 




(3.1) 




(3.2) 



The spatial curvature of the LTB geometry is 




(3.3) 
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For our choice of f(r) we find that at the initial time 

^^-^t 1 -^)- 6 ^! 1 -^)' (3 - 4) 

Overdense regions have positive spatial curvature, while underdense ones negative 
curvature. This is very similar to the initial condition considered in the model of spherical 
collapse [32]. 

When the inhomogeneity is denser near the center, we have f(r) < for r < ro and 
f(r) = for r > ro. It is then clear from eq. (2.4) that, in an expanding Universe, the 
central region will have R = at some point in its evolution and will stop expanding. 
Subsequently, it will reverse its motion and start collapsing. The opposite happens if the 
inhomogeneity has a central underdensity. In this case, the central region expands faster 
than the surrounding denser spherical shell. The width of the shell decreases, while its 
density increases. It is the latter configuration that is relevant if we want to model the 
Universe as being composed mainly of voids separated by thin dense regions. 

Our expressions simplify if we switch to dimensionless variables. We define i = tHi, 
f = r /r , R = R/r , where Hf = p i /(QM 2 ) is the initial homogeneous expansion rate and 
ro gives the size of the inhomogeneity in comoving coordinates. The evolution equation 
becomes 

R 3M(f) /(f) 

R? ~ AnR 3 + R? ' ^ 
with M = M/(p 0)i rl) and / = 6M 2 f /(p 0:i r%) = f/Hf, Hi = H^. The dot now denotes 
a derivative with respect to t. 

We take the affine parameter A to have the dimension of time and we define the 
dimensionless variables A = HiX, k° = k°, k 1 = k 1 /Hi, k 3 = rok 3 . The geodesic equations 
(2.6)-(2.8) maintain their form, with the various quantities replaced by barred ones, and 
the combination 1 + / replaced by H^ 2 + /. For geodesies going through subhorizon 
perturbations with ffj < 1 the effective curvature term / plays a minor role. However, 
this term is always important for the evolution of the perturbations, as can be seen from 
eq. (3.5). The optical equation takes the form 

TT^r = ~A k )> (3 - 6) 

with p = p/po,i- We omitted the shear, as it gives a negligible contribution to our results. 
The initial conditions (2.14), (2.15) become 

d\[A 

H 



dX 
A 



A=0 

with A = H 2 A and Cl = H 2 VL. 



V s = Jn s . (3.7) 
0, (3.f 
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Figure 1. The evolution of the density profile for a central underdensity surrounded by 
an overdensity. 



4. Single crossing 



The typical cosmological evolution is displayed in fig. 1 for a central underdensity that 
is surrounded by an overdense region. The initial density pi(r) = Pi(r) / p 0)i = 1 + e(r) is 
constant pi = 1 + ei in the region f < 0.8, constant pi = l + €2 in the region 0.9 < f < 0.95, 
and pi — 1 for f > 1. In the intervals 0.8 < f < 0.9 and 0.95 < f < 1 it interpolates 
linearly between the values at the boundaries. This guarantees that the conditions 
f(l) = f'(l) = 0, imposed by the matching of the FRW and LTB metrics, are satisfied. 
We use €\ = —0.01. The value of €2 is fixed by the requirement that Jq 1 e(f)f 2 df = 0. In a 
previous publication [27], we studied in detail the evolution of such inhomogeneities. We 
showed analytically that their growth is consistent with the standard theory of structure 
formation and the spherical collapse model [32]. Inhomogeneities that are initially of 
horizon size (Hi ~ 1) with |ei| ~ 10~ 5 and have a characteristic scale O(10) h~ l Mpc 
today also have a density contrast of 0(1). In our numerical solutions we do not follow 
the very early evolution of the inhomogeneities. The perturbations we consider are already 
of subhorizon size and have a density contrast of O(10~~ 2 ). They evolve to form present- 
day structures with size O(10) h" 1 Mpc or larger and density contrast of 0(1). The 
profile of the perturbations that we assume in this work differs slightly from the one in 
ref. [27] with respect to the width of the shell. In ref. [27] the widths of the underdense 
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Figure 2. The relative difference in redshift, coordinate time, and beam area, between the 
propagation in the background of fig. 1 and in a homogeneous background, as a function 
of the radial coordinate. 



and overdense regions were taken equal. In this work we assume that the initial shells 
are narrower and denser than in ref. [27]. Clearly, our modelling of the Universe cannot 
reproduce all the details of the statistical nature of the inhomogeneities, especially during 
their evolution in the non-linear regime. Our choice results in a model of the Universe 
dominated by voids, in rough agreement with observations. It is possible, however, that 
the effect of underdensities is overestimated, as they occupy a slightly larger fraction of 
the total volume than the overdense regions already in the initial small perturbations. 

In fig. 1 we display the density profile at times £=10, 100, 200, 250. We follow the 
evolution at later times as well, even though we do not depict it in fig. 1. We normalize the 
energy density to that of a homogenous FRW background (given by f>FRw(t) = p(t, 1))- 
The earliest time corresponds to the curve with the smallest deviation from 1, while 
the latest to the curve with the largest deviation. We observe that the density contrast 
grows and eventually becomes of 0(1). The central energy density drops relative to the 
homogeneous background, while the surrounding region becomes denser. The radius of the 
central underdensity grows relative to the total size of the inhomogeneity, as this region 
expands faster than the average. The surrounding shell becomes thinner and denser. 
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The effect of the spherical inhomogeneity on the characteristics of a light beam is 
depicted in figs. 2, 3. (Fig. 3 is a magnification of fig. 2 around f = 0.) We assume a 
perturbation with Hi = 1 at the initial time ij = 0, so that the effects on the light beam 
are clearly visible. We use this perturbation for the numerical analysis both in this section 
and in section 5. Such a perturbation leads to an inhomogeneity at the present time with 
size that would be in conflict with observations (approximately 350 h -1 Mpc). We consider 
realistic perturbations in section 6. We consider a beam with = that is emitted from 
a point with f s = 1.5 and passes through the center of symmetry. The emission time is 
i s = 250 and the signal reaches f = 1.5 again at a time i ~ 441. The redshift for the exiting 
beam at f = 1.5 is z ~ 0.46. We plot the relative difference in redshift (z — zfrw) I 'zfrw, 
coordinate time (i—iFRw)/tFRW, and beam area (^/A — \J A FRW ^j / \j A FRW , between the 
propagation in the background of fig. 1 and in a homogeneous background, as a function 
of the radial coordinate f. The arrows indicate the evolution of various quantities as the 
beam enters the inhomogeneity from one side and exits from the other. In the region 
0.1 ^ f ^ 0.9 we observe a significant deviation of all quantities from their values in 
a homogeneous background. However, the deviation becomes small in the region near 
the origin. When the beam moves out of the inhomogeneous region only the beam area 
deviates from the value in a homogeneous background. The coordinate time and the 
redshift are not affected significantly. 

We can obtain an understanding of the evolution depicted in figs. 2, 3 through an 
analytical treatment. The effect of the inhomogeneity on the characteristics of the light 
beam can be estimated analytically for perturbations with size much smaller than the 
distance to the horizon. These have Hi <C 1. Let us consider a beam with = that 
passes through the center of the spherical inhomogeneity. The null condition (2.7) can be 
written as 



neglect the second term in the above expression. This indicates that the spatial curvature 
does not play a role if the inhomogeneities are much smaller than the horizon. We can 



also employ the approximation R'(i,r) ~ R'(i s ,f) + R (i s ,r)(i — t s ), as the time it takes 
for the light to cross the inhomogeneity is much shorter than the Hubble time. In fact, 
t — t s = O(Hi) (see eqs. (4.2), (4.3) below). We denote by f s the location of the source 
and by t s the emission time of the beam. The solution of eq. (4.1) is 





for incoming and outgoing beams respectively. If we keep terms up to O(Hf) we can 



i-i s = Hi (R(t 8 , f s ) - R{t„, f)) + Hi 



2 



R'(t s ,r)R(t S: f)df 




(4.2) 
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Figure 3. The relative difference in redshift, coordinate time, and beam area, between the 
propagation in the background of fig. 1 and in a homogeneous background, as a function 
of the radial coordinate. 

for an incoming beam, and 

t-t„ = Hi (R(i s , f) - R(t„ r 8 j) - Hf [ R'(i s , f)R(t s , f)alf 

+ Hf (R(t s , f) - R(t„ f s )) R(t„ f) + 0(H?) (4.3) 

for an outgoing one. These expressions are confirmed by numerical solutions. 

We can also make a comparison with the propagation of light in a FRW background. 
In this case we have R(t, f) = a(t)f = R(i, l)f. We have expressed the scale factor in 
terms of the value of the function R(t, f) at the boundary of the inhomogeneous region 
f = 1. Let us consider light signals emitted at f s — 1 and observed at the center (f G = 0) 
of the inhomogeneity. The difference in propagation time within the LTB and FRW 
backgrounds is 

to-(t ) FRW = Hf £ft(t.,r)R(t e ,r)dr-^R(i^^ 

For signals originating at f s — and detected at f Q = 1 the time difference has the 
opposite sign. As a result, the time difference for signals that cross the inhomogeneity is 
of 0(Hf). This is in agreement with figs. 2, 3. 
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We can derive similar expressions for the redshift of a light beam that passes through 
the center of the inhomogeneity. The geodesic equation (2.6) can be written as 



for incoming and outgoing beams respectively. In this way we find 
ln(l + z)=H i {Rifs, f.) - R(t„ f)) 

- f R{t„ f) (R(t s , f s ) - R(i SJ fj) air + O(Hf) (4.6) 

Jr s v 7 

for an incoming beam, and 

m(l + z)=Hi (R(t s , f) - R(i s , f s )) 

+ Hi f R(i s , f) (R(t„ f) - R(i s , f a )) air + 0(H?) (4.7) 

Jr B y ' 

for an outgoing one. These expressions are confirmed by numerical solutions. 
For signals originating at f s — 1 and detected at f — the redshifts obey 

In ( ] + Z ) = Hf f 1 R(t„ f) (R(t B , 1) - R(t s , f)) df 
\l + z FRW J Jo v 7 

- ^R{t s , l)R(t 8 , 1) + O(Hf). (4.8) 

For signals originating at f s = and detected at f = 1 the r.h.s. of the above equation 
has the opposite sign. As a result, the redshift difference for signals that cross the 
inhomogeneity is of O(Hf). Again, this is in agreement with figs. 2, 3. 

5. Luminosity distance and flux conservation 

The beam area obeys the second-order differential equation (3.6), whose solution depends 
crucially on the initial conditions. For initial conditions given by eqs. (3.7), (3.8) and 
symmetric situations, it is possible to determine the solution analytically. For signals 
emitted from some point f s at a time t = t s and observed at f — we have [33, 16] 

VI= (1 + z)R(i s ,f s )Vfl. (5.1) 

This is in agreement with figs. 2, 3. For signals emitted from the center f s = and 
observed at f a at a time i Q we have 

\ / A = R(i ,f )Vn. (5.2) 

However, for a signal that crosses the inhomogeneity we need to integrate eq. (3.6) from 
f = to f a with initial conditions determined by the propagation from f s to f — 0. These 
are different from (3.7), (3.8), so that an analytical solution is not easy. 

Using perturbation theory, it is possible to obtain an analytical estimate of the 
deviation of the luminosity distance from its value in homogeneous cosmology. We 
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Figure 4. The deviation of the luminosity distance from its value in a homogeneous 
background, as a function of the impact parameter. 



consider beam trajectories that start at the boundary of the inhomogeneity, pass through 
its center and exit from the other side. The optical equation (3.6) can be written in the 
form 



1 -^^+ - . =--p k°) VA. 




k )_^ + lx^F = -2 ii ( k ) VA - (5 ' 3) 

We express dki/dX in the above equation using the geodesic equation and obtain 

2HM_ _1T Hff \ dVA 3 r 
-Jjff R> + 2(1 + Hff)) df ~ 2 P l + HU [b - } 

where the positive sign in the second term corresponds to ingoing and the negative sign 
to outgoing geodesies. 

We use a simplified initial configuration for this estimate. We take p(0, f) = for 
f <f\ and p(0, f) = 1/(1 — f\) for f > f\. The initial conditions for an ingoing beam can 
be taken vl(f = 1) = 0, dvA(f = l)/df = —1, without loss of generality. We use the 
expansion 

y/A = VW + HiVAW + H?VAW + 0(H?), (5.5) 



and calculate v AW in each order of perturbation theory. The calculation in described 
in the appendix. We point out that our choice of initial configuration, that involves a 
discontinuous energy density, results in the appearance of 5-function singularities in the 
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second derivatives of R with respect to f. These must be taken into account in a consistent 
calculation, as described in the appendix. We have checked that the expressions (5.1) and 
(5.2) are reproduced correctly by our results, up to second order in Hi. 
When the photon exits the inhomogeneity at f = 1 we find 

VA(f=\) = 2 + 4i?, + (5 - - 3 - \ H\ + O(Hf) (5.6) 



= !) = ! + m, + ( 6 - Hf + Cm (5.7) 

For a homogeneous universe we have 

y/A(r = 1) = 2 + 4i^ + 2Hf (5.8) 

^?(r = 1) = l + 4jZ i + 3ff?. (5.9) 
ar 

It is clear that the effect of the inhomogeneity is of O(Hf). This conclusion has 
been confirmed by numerical solutions of the exact optical equations, without any 
approximations. Several beam crossings were studied for a multitude of values of Hi. 
The coordinate time t needed for the crossing, the redshift z and the beam area A were 
plotted in terms of Hi. The polynomial fits to these plots verify with good precision 
the conclusions of the present and previous sections: The deviations of t Q and z from 
their values in a homogeneous background are of O(Hf), while the deviations of A are of 
0{H1). 

From the above we can conclude that the main effect of an inhomogeneity is to 
modify the luminosity distance of a light source behind it, while leaving the redshift 
largely unaffected. In fig. 4 we depict the modification of the luminosity distance 
{D L - (D L ) frw)/(Dl) frw for light beams that cross a spherical inhomogeneity at 
various angles. The inhomogeneity is the one employed in the previous section. Its size 
is approximately 350 /i -1 Mpc, much larger that what is deduced from observations. We 
consider realistic inhomogeneities in the following section. The crossing at a varying angle 
is achieved by choosing different values for the constant in eq. (2.8). The light source 
is always located at the same point with f s = 1.5. The beam is allowed to propagate until 
the redshift reaches a certain value z Q ~ 0.46. The beam with = that goes through 
the center reaches r a ~ 1.5 at this time. The difference of the resulting luminosity distance 
from the one in a FRW background is depicted in fig. 4 as a function of = c^jc^^ max . 
The largest value c^ maa; corresponds to beams that are emitted tangentially with respect 
to the center of symmetry. The angle 9 between the initial direction of the beam and the 
radial direction is determined through the relation sin# = c^. The variable plays the 
role of an impact parameter, normalized to 1 for light beams emitted tangentially. In fig. 
4 we observe that, if is sufficiently small for the light to travel through the central 
underdense region, the luminosity distance is enhanced relative to the homogeneous 
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case. For larger values of the light travels mainly through the overdense shell and 
the luminosity distance in reduced. 

If the redshift is not affected significantly by the propagation in the inhomogeneous 
background, the conservation of the total flux implies that the average luminosity distance 
must be the same as in the homogeneous case [34, 35]. As we have seen, this is not the 
case for an observer located at the center of an underdensity. The resulting increase in the 
luminosity distance, arising mainly from the increase in the redshift, can by employed for 
the explanation of the supernova data, even though the size of the required inhomogeneity 
is probably in conflict with the observed large-scale structure [15]-[24]. On the other hand, 
if the source and the observers are located outside the inhomogeneity, as for fig. 4, so that 
the redshift is essentially unaffected, we expect that the energy flux may be redistributed 
in various directions but the total flux will be the same as in the homogeneous case [34, 35]. 
This implies that the integral of (Dl — (Dl) frw ) / (Dl) frw over all angles must vanish. 
For an isotropic source the integration over the solid angle is 2tt sin 9 d9. This is equivalent 
to the integration over the impact factor c^, with a weight c^/Jl — c|. The integral of 
the function depicted in fig. 4 is indeed approximately zero. An analytical proof is very 
difficult, but a numerical analysis shows that a cancellation at the 90% level takes place 
between the positive and negative contributions. The remaining deviation from zero is 
caused by numerical errors and the small (of O(Hf)), but non- vanishing, difference in 
redshifts in inhomogeneous and homogeneous backgrounds. 

6. Multiple crossings 

In this section we consider light beams that pass through several of the inhomogeneities 
described in the previous sections. The light is emitted at some time i s from a point with 
f = 1 at the edge of the inhomogeneous region. Its initial direction is assumed to be 
random. The initial conditions for the beam area are given by eqs. (2.14), (2.15). The 
light moves through the inhomogeneity and exits from a point with f = 1. Subsequently, 
the beam crosses the following inhomogeneity in a similar fashion. The angle of entry 
into the new inhomogeneity is assumed to be random again. The initial conditions are 
set by the values of \f~A and d\fA/d\ at the end of the first crossing. This process is 
repeated until the light arrives at the observer. Of course, as time passes the profile of 
the inhomogeneities changes, as depicted in fig. 1. 

The assumption of entry into a spherical inhomogeneity at a random azimuthal angle 
is realized by selecting the impact factor with a probability ~ c^dc^. The absence of 
the denominator ^Jl — c| that appears in the weight employed in the previous section is 
justified by the fact that the source is located far from the center of the inhomogeneity, so 
that c^Cl. This is not strictly true for the first 1 or 2 crossings, but the induced error 
is small. In order to place the observer outside the inhomogeneities we always replace the 
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Figure 5. The distribution of luminosity distances for various redshifts in the LTB 
Swiss-cheese model if the inhomogeneities have a characteristic scale o/40/i _1 Mpc. 



final spherical region crossed by the beam with a homogeneous configuration. 

The essence of our procedure is that the light beam encounters various 
inhomogeneities at various angles along its path. In this section we consider the 
propagation of the beam only in the intervals < r < 1 of coordinate systems with origins 
at the center of the various inhomogeneities. This means that we neglect the homogeneous 
region between the inhomogeneities that we assumed in the previous region. (Essentially 
we assume that its width is negligible.) The reason for this omission is that the presence 
of the homogeneous region generates a bias towards negligible deviations of the luminosity 
distance from its value in a homogeneous cosmology. For example, if we consider light 
propagation in the intervals < f < 1.5 as in the previous section, a large number of 
beam trajectories propagate only within the homogeneous regions with 1 < f < 1.5. 
These give a luminosity distance for the source equal to that in the FRW cosmology. On 
the other hand, there are trajectories that propagate through the inhomogeneities, for 
which the cancellation between positive and negative contributions results in a negligible 
total deviation of the luminosity distance from the value in a homogeneous background. 
It is the latter events that we are interested in, while the former are rather unphysical. 
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Figure 6. Same as in fig. 5 for a characteristic scale of 133 h 1 Mpc. 



The procedure we outlined above has an unsatisfactory element. The spherical 
inhomogeneities are assumed to follow each other continuously along the beam trajectory. 
If the angles of entry of the beam are non-zero, the resulting geometry implies that there 
is an overlap of the inhomogeneities outside the beam trajectory. This problem cannot 
be corrected as long as the assumption of spherical symmetry of the inhomogeneities is 
maintained. A choice must be made between an artificial bias in the luminosity distance if 
intermediate homogeneous regions are introduced, or the overlap of the inhomogeneities 
outside the beam trajectory. We have chosen the second option, as we believe that it 
provides a more reliable estimate of the distribution of luminosity distances. 

The total number of crossings determines the redshift and the final beam area, related 
to the luminosity distance. We determine the arrival time at the observer by requiring that 
the redshift reach a specific value. We repeat the calculation many times (at least 1000) 
for each value of the redshift and plot the resulting distribution of luminosity distances. 
The emission time is such that the arrival time is t ~ 284 for all the redshifts that we 
consider. At this time the profile of the inhomogeneity is very similar to the curve in fig. 
1 with the largest deviation from 1. The deviations of the exact arrival time t Q from the 
time in a homogeneous background is one order of magnitude smaller than the respective 
deviation for the luminosity distance. This is in agreement with the discussion in the 
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Figure 7. Same as in fig. 5 for a characteristic scale of 400 h 1 Mpc. 



previous section. 

In figs. 5, 6, 7 we depict the distributions of the deviations of the luminosity distances 
from the value in a homogeneous background for various redshifts. The three figures 
correspond to inhomogeneities with different characteristic length scales at the time of 
the emission of light. The background through which the light propagates is constructed 
as described in the previous sections. At the initial time tj = the inhomogeneities have 
the profile descibed in the previous section with e\ = —0.01. The subsequent evolution 
is depicted in fig. 1. The characteristic length scale of the inhomogeneities relative to 
the distance to the horizon at the time ti — is ro/(i?, rl ) = Hi. This quantity does 
not appear in the rescaled evolution equations for the background. It appears only in 
the rescaled geodesic equations. For this reason we can use the same background, with 
an evolution depicted in fig. 1, in order to discuss inhomogeneities of various length 
scales. The important phenomenological quantity is the scale of the inhomogeneities 
today. This is given by R(tf, r )/(Hj x ) = R(tf,r ) = R(if,l)Hi. The rescaled present 
time if is equal to the time of arrival of light signals to the observer i Q for all the cases 
we consider. As we mentioned already, if — i a ~ 284. Our solution has R(i ) = 0.133. 
Using Hj 1 = 3 x 10 3 h~ l Mpc we have R(t f , r ) = 400 h~ l Hi Mpc. 
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The three figures 5, 6, 7 correspond to Hi = 1/10, 1/3, 1, respectively. They describe 
the effect on the propagation of light of inhomogeneities with sizes 40, 133, 400 h^ 1 Mpc 
today. The present profile of the inhomogeneities has a density contrast 0(1). Their 
evolution, as modelled by the LTB metric, is roughly consistent with the standard theory 
of structure growth. The choices Hi = 1/3 and 1 lead to present-day inhomogeneities with 
length scales larger than those in typical observations. The size of the same perturbations 
at horizon crossing is larger than the value ~ 10 -5 implied by the CMB. However, we have 
included them for two reasons. Firstly, because there are indications that the presence 
of such large structures may be supported by observations [25, 26]. Secondly, because we 
would like to understand if inhomogeneities with sizes comparable to the horizon distance 
can have a significant effect on the luminosity distance for a random location of the 
observer. 

The total integral of the distributions has been normalized to 1 in all cases, so that 
they are in fact probability densities. They have similar profiles that are asymmetric 
around zero. Each distribution has a maximum at a value larger than zero and a long 
tail towards negative values. The average deviation is zero to a good approximation in all 
cases. This is expected according to our discussion of flux conservation in the previous 
section. As long as the light propagation in an inhomogeneous background does not 
modify significantly the redshift, the energy may be redistributed in various directions, 
but the total flux is conserved and remains the same as in a FRW background. 

The longer tail of the distribution towards luminosity distances smaller than the 
one in a homogeneous background is a consequence of the presence of a thin and dense 
spherical shell around each central underdensity. The number of beam trajectories that 
propagate through several shells is small. However, the focusing of the beam is substantial 
for such beams and the resulting luminosity distance much shorter than the average. The 
effect of the long tail is compensated by the shift of the maximum of the distribution 
towards positive values. The form of the distribution is very similar to that derived in 
studies modelling the inhomogeneities through the standard Swiss-cheese model [11]. In 
that case the strong focusing is generated by the very dense concentration of matter at 
the center of each spherical inhomogeneity. We emphasize, however, that the two models 
have a different region of applicability. The standard Swiss-cheese model is appropriate 
for length scales of 0(1) h~ l Mpc or smaller, while the LTB Swiss-cheese model for scales 
of 0(10) h~ x Mpc or larger. 

7. Determination of cosmological parameters 

The form of the distributions can be quantified in terms of two parameters: The width of 
the distribution 5d and the location of its maximum 8 m > 0. The first one characterizes 
the error induced to cosmological parameters derived through the curve of the luminosity 
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Figure 8. The relative deviation of the luminosity distance from its value for w = 
(non-relativistic matter) in homogeneous cosmology. 

distance as a function of redshift, while the second one the bias in such determinations. 
As discussed extensively in ref. [11], a small sample of data is expected to favour values of 
the luminosity distance near the maximum of the distribution, and thus generate a bias. 
In figs. 5, 6, 7 we observe that both Sd and S m grow with inreasing redshift z and scale 
Hi. For Hi = 1/10 the average width Sd increases from approximately 0.005 to 0.01 as 
z increases from 0.5 to 2. The maximum is 8 m ^ 0.002 for all z. The asymmetry of the 
distribution is very small. For Hi = 1/3 the average width increases from approximately 
0.005 to 0.02 as z increases from 0.5 to 2. The maximum is 5 m ^ 0.002 for all z. The 
asymmetry of the distribution and the longer tail towards negative values are clearly 
visible in fig. 6. For Hi = 1 the average width increases from approximately 0.01 to 0.04 
as z increases from 0.5 to 2. The maximum increases from 0.005 to 0.01. The asymmetry 
of the distribution is very distinctive in fig. 7. 

The values of 5d and 5 m are very small for all the values of Hi that we considered. 
This implies that we do not expect a significant effect on the cosmological parameters. 
As an interesting example we consider the parameter w that appears in the equation of 
state of the cosmological fluid. In homogeneous cosmology the luminosity distance is a 
function of w and the redshift z. The relative deviation of the luminosity distance from 




Figure 9. The function a{z) = — ds(w,z)/dw\ 



its value for w = (non-relativistic matter) is 

( \ d l{w,z) i i - + ir (1+3w)/2 1 ryn 

S (W, Z) = — — : 1 = 1. 7.1 

1 ' D L (0,z) l + 3w i_( z + i)-V2 ^ > 

We depict this function in fig. 8. For w — — 1 the luminosity distance is larger by roughly 
35% relative to w = for 2 = 0.5. For z = 1 the relative increase is approximately 70%. 

The deviations of 5d(z) and 6 m (z) from zero because of the presence of 
inhomogeneities can be attributed to deviations of w from zero if the cosmology is assumed 
to be homogeneous. This can be achieved by identifying Sd(z) or S m (z) with s(w,z). For 
small w we have 

s(w, z) ~ -3 [ 1 - It—^TT-t) w = -<x(z)w- (7-2) 



;z + i)V2_i ; 

We depict the function a(z) in fig. 9. For a given value of z we can derive effective values 
of w from the relations w e ff = —5d, m (z)/a(z). It is clear from the values we quoted above 
for 5d, S m and fig. 9 that |iw e //| C 1. As a result, for a random location of the observer 
the perceived acceleration of the Universe cannot be attributed to the modification of the 
luminosity distance by large-scale inhomogeneities, even when their characteristic scale is 
smaller than the horizon distance by less than a factor of 10. 

On the other hand, the presence of inhomogeneities induces a statistical error in 
the value of w deduced from astrophysical data, as well as a shift of its average value 
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if the sample is small. According to our results, for Hi = 1/10 (and present size of the 
inhomogeneities 40 hr x Mpc) the error is 5w ~ 0.015 for all z between 0.5 and 2, while 
the average value w for a small sample is negative and of O(10~ 3 ). For Hi = 1/3 (and 
present size of the inhomogeneities 133 h -1 Mpc) the error increases from 0.015 to 0.025 
as z increases from 0.5 to 2, while the average value w is again negative and of O(10~ 3 ). 
For Hi = 1 (and present size of the inhomogeneities 400 h~ x Mpc) the error increases from 
0.03 to 0.05, while the average is w ~ —0.015. The shift in the average value is always 
smaller than a standard deviation. 

The values of 5d and 5 m can be compared to those generated by the effects of 
gravitational lensing at scales typical of galaxies or clusters of galaxies. At such scales 
the Universe is modelled through the standard Swiss-cheese model, with the mass of each 
inhomogeneity concentrated in a very dense object at its center [11]. The typical values 
of 5d and 5 m are larger by at least an order of magnitude than the ones we obtained. The 
reason is that in our model the density contrast is always of 0(1). Our study indicates that 
the effect on the luminosity distance grows when the length scale of the inhomogeneities 
increases and becomes comparable to the horizon distance. However, the strong lensing 
effect of a high concentration of mass, such as a galaxy or cluster of galaxies, gives a much 
larger effect. 

We conclude that within our model the presence of inhomogeneities with large length 
scales, even comparable to the horizon distance, and density contrast of 0(1) does not 
influence significantly the propagation of light if the source and the observer have random 
locations. It is possible that our modelling of the Universe is lacking some essential feature 
that could generate a significant effect on the luminosity distance. For example, it has 
been suggested that large fluctuations in the spatial curvature may result in a strong 
backreaction on the overall expansion [37]. Unfortunately, there are no known exact 
models that realize such a scenario. In their absence the effect on the luminosity distance 
cannot be computed. In our model, there are variations of the local curvature which result 
in the collapse of the overdense regions (either central overdensities, or shells surrounding 
underdensities). The evolution is in approximate agreement with the standard theory 
of structure formation. It seems difficult to reconcile a much larger local curvature 
with a large-scale structure consistent with observations. Another possibility is that the 
assumed spherical symmetry of the inhomogeneities in our modelling of the Universe is 
too constraining. It is possible that photon propagation through inhomogeneities without 
such a symmetry results in a stronger modification of the luminosity distance. This point 
will be the subject of a future investigation. 
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8. Appendix 

In this appendix we describe the solution of eq. (5.4) through the use of perturbation 
theory in Hi. In this way we demonstrate that a spherical inhomogeneity induces a 
deviation of the luminosity distance from its value in a homogeneous background which 
is of 0(Hf). 

The null constraint (2.7), keeping terms up to O(Hi), is di/df = ±Hi, with the 
negative sign corresponding to ingoing and the positive to outgoing geodesies. We can 
set t s = so the geodesic inside the inhomogeneity is 

t = -Hi(r-r ) + O(H?) (8.3) 
for ingoing, and 

i=Hi{r + f ) + O(H?) (8.4) 

for outgoing geodesies. We can treat t as an O(Hi) quantity. 

The initial configuration we use for this estimate has p(0, f) = for f < f\ and 
p(0,f) = 1/(1 — ff) for f > f\. From (3.5) we can calculate various derivatives of R at 
i=0: 

R(t, f) = R(0, f) + tR(0, f) + 0(H 2 ), (8.5) 

§(t,f) = ^R(0,f) + O(Hf). (8.6) 

For f > f i we have 

R(0,f) = + 1 (8.7) 
2 r i ^ _ xj 



/2 (0,f) 



3 



3ff 



r 4 (ff - 1 



For r < T\ both i? (0, r) and i? (0, r) are zero. For the initial configuration that we 
assume, R is a continuous function of f. However, R is discontinuous at f = fx and 
f = 1, while i? has 5-function singularities at the same points. 

The initial conditions for the solution of eq. (5.4) for an ingoing beam can be taken 
\A(1) = 0, dy/~A(l)/df = — 1, without loss of generality. We use the expansion 

VX=VM + HiVW+i%VM + 0(H?). (8.9) 

To zeroth order in Hi, eq. (5.4) becomes d 2 V A^/df 2 = 0, with solution V AW (f) = 
— (r — 1) for ingoing and V AW (f) = r + 1 for outgoing beams. To first order in Hi, 
eq. (5.4) gives dy/A^/df = —2, with solution V A^ l \f) = r 2 — 2r + 1 for ingoing and 
vIW (f) = r 2 + 2r + 1 for outgoing beams. 



Light Propagation and Large-Scale Inhomogeneities 



25 



To second order in Hi, and for r > r\ and ingoing geodesies, we obtain 



= - -p(0,f)VAW. (8.10) 

For f <f\ and ingoing geodesies, we have 

ol 2 VW) + f{f)dVW +2 c/v / I« _ 
(if 2 2 df df 

For f < fi and outgoing geodesies, we have 

d 2 ^W) f'(r)dVAW ^ d^W) 
df 2 2 df df 

Finally, for f > f\ and outgoing geodesies, we obtain 



___ + / _2t( f )i2 (0, f) - -i? (0, f) + — — - 2 



= - -p(0,f)Vl(o). (8.13) 

The above equations can be solved analytically through simple integration, with the 
values at the end of each interval determining the initial conditions for the next one. 
The only non-trivial point is that the 5-function singularities of R at f — f i and f = 1 
induce discontinuities in the values of d\J AW) / df at these points. These must be taken 
into account in a consistent calculation. The discontinuities can be easily determined 
through the integration of eqs. (8.10)-(8.13) in an infinitesimal interval around each of 
these points. The remaining calculation is straightforward. When the photon exits the 
inhomogeneity at f = 1 we find that \/~A is given by eq. (5.6). 
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